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We present a model which describes the properties of odd-even nuclei with one nucleon more, or 
less, with respect to the magic number. In addition to the effects related to the unpaired nucleon, 
we consider those produced by the excitation of the closed shell core. By using a single particle basis 
generated with Hartree-Fock calculations, we describe the polarization of the doubly magic-core with 
Random Phase Approximation collective wave functions. In every step of the calculation, and for all 
the nuclei considered, we use the same finite-range nucleon-nucleon interaction. We apply our model 
to the evaluation of electric quadrupole and magnetic dipole moments of odd-even nuclei around 
oxygen, calcium, zirconium, tin and lead isotopes. Our Random Phase Approximation description 
of the polarization of the core improves the agreement with experimental data with respect to 
the predictions of the independent particle model. We compare our results with those obtained 
in first-order perturbation theory, with those produced by Hartree-Fock-Bogolioubov calculations 
and with those generated within the Landau-Migdal theory of finite Fermi systems. The results of 
our universal, self-consistent, and parameter free approach have the same quality of those obtained 
with phenomenological approaches where the various terms of the nucleon-nucleon interaction are 
adapted to reproduce some specific experimental data. A critical discussion on the validity of the 
model is presented. 

PACS numbers: 21.10.Ky ; 21.60.Jz 


I. INTRODUCTION 

The description of the angular momenta of odd-even nuclei is one of the greatest successes of the nuclear shell-model. 
The angular momentum and the parity of the single particle (s.p.) level of the unpaired nucleon correspond to those 
of the whole nucleus. This extreme shell model picture, where the s.p. properties are imposed to the whole interacting 
many-body system, is weakened when observables which are not quantized are investigated. For example, the values 
of electric quadrupole and magnetic dipole moments are rather different from those predicted by the extreme shell 
model. For these quantities, the interaction between the unpaired nucleon and the other nucleons plays an important 
role, strongly modifying the pure shell model predictions. 

One expects that in odd-even nuclei around doubly-magic ones, the effects related to the extreme shell model 
picture, which we shall name henceforth Independent Particle Model (IPM), are dominant with respect to those 
induced by the two-body part of the nuclear hamiltonian, usually called residual interaction. This has induced the 
development of various perturbative models [1, 2]. A class of these models is based on the straightforward application 
of the traditional perturbation theory, which, however, has convergence problems, at least for the terms beyond those 
of the first order [3]. 

Another perturbative approach consists in considering the residual interaction within the framework of the Random 
Phase Approximation (RPA) theory. This approach is based on the extension of the Landau theory of quantum liquids 
to finite Fermi systems (FFS) done by Migdal [4]. This theory has been applied with great success to describe the 
properties of odd-even nuclei in the region of ^*^®Pb [5-7]. In these works the s.p. wave functions are generated 
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by diagonalizing a Woods-Saxon well, and the residual interaction is a zero-range Landau-Migdal force [4]. More 
recently, this approach has been extended to use the same interaction in the production of the s.p. wave functions and 
in the RPA calculations [8-12]. Furthermore, also pairing effects have been considered. This self-consistent theory 
of FFS has been applied also to the tin and to the calcium isotopes. In these approaches, the use of the concept of 
quasi-particles implies the definition of the effective charge which contains free parameters whose values are chosen 
to have a good description of the data. In the phenomenological approach of Refs. [5-7] also the constants defining 
the strength of the various channels of the Landau-Migdal nucleon-nucleon effective interaction have been used as 
free parameters. In this case, the values of these parameters are related to the dimension of the s.p. configuration 
space. The application of these methods to a different region of the nuclear chart, or, more simply, to a different 
s.p. configuration space, requires a new selection of the free parameter values. The self-consistent FFS approach of 
Refs. [8-12] does not have cut-off problems in the RPA calculations, since it considers the full s.p. configuration 
space, continuum included. The problems arise in the pairing calculations, where the free parameter is the value of 
the cutoff energy used as regulator, in analogy to the renormalisation procedure adopted in effective field theories. 

Inspired by these works, we propose here an extension of the RPA formalism to construct a fully self-consistent 
approach. For each nucleus considered, we use a set of s.p. wave functions generated by a Hartree-Fock (HF) 
calculation. The same, effective, nucleon-nucleon interaction used in HF is adopted in the RPA calculations, which 
allows the evaluation of the odd-even nucleus observables. The parameters of the force, which we consider having 
finite-range, have been chosen once forever in a global HF fit of properties mainly related to the ground state of a large 
set of nuclei in all the regions of the nuclear chart [13]. The use of a finite-range interaction ensures the stability of our 
results whose values, after convergence has been reached, do not depend any more on the size of the s.p. configuration 
space. 

The basic ideas of our work are relatively simple and straightforward. First, we use effective nucleon-nucleon 
interactions whose parameters have been chosen to reproduce data which are not directly related to the observables 
we investigate. Second, the chain of calculations needed to arrive at the final results does not contain any additional 
free parameter, but they depend exclusively on the effective nucleon-nucleon interaction. To be sure of identifying 
features related to the theory, and not to the specific choice of the interaction or to the nucleus investigated, we 
carried out calculations with different nucleon-nucleon forces, for a set of nuclei in a large region of the nuclear chart. 
We mainly concentrate our investigation on two observables, the electric quadrupole, Q, and the magnetic dipole, /r, 
moments of a selected set of odd-even nuclei. The evaluation of these quantities requires, respectively, the description 
of the 2+ and 1+ excitations of the even-even core nuclei which we obtain by using a RPA approach [14]. 

In Sec. H we present the basic ideas our model. In this section we also specify the form of the electromagnetic 
operators that we have adopted to describe Q and fj,. The expression of the operator for Q is simple and straight¬ 
forward. More complex is the situation for /r, where, in addition to the traditional one-body operator, we consider 
also two-body currents generated by the exchange of charged pions. We shall refer to these latter currents as meson 
exchange currents (MFC). In Sec. HI we discuss some details of the calculations related to the specific applications of 
the model. The results obtained for Q and /r are presented in Sect. IV. In both cases we also study the excitation of 
the even-even core nuclei for the two multipolarities involved. In Sec. V we summarize the main results of our study, 
and we draw our conclusions. 


II. THE MODEL 

The starting point of our model is the construction of the basis of s.p. states \(l)a), which we generate by solving 
the HF equations. We used the symbol a to indicate all the quantum numbers characterizing the s.p. state, i.e. 
the principal quantum number Uq, the orbital angular momentum Iq,, the total angular momentum j^, its third 
component nia, and the third component of the isospin ta = —1/2 for protons, and ta = 1/2 for neutrons. 

In the IPM picture, all the s.p. states below the Fermi energy are completely occupied, and those above it are 
empty. We indicate with |$o) the Slater determinant describing the IPM ground state of the doubly magic nucleus 
composed by A nucleons. By definition, the RPA ground state of the doubly magic nucleus, which we indicate as 
|4'o), contains correlations beyond the IPM. 

The odd-even nuclei we want to describe are obtained by adding, or subtracting, one nucleon to the doubly-magic 
A-nucleon system. We label the states of these nuclei as | A ± I; a), where the symbol a indicates the set of quantum 
numbers characterizing the s.p. state of the added, or subtracted, nucleon. 

We describe the response of the odd-even nucleus to the perturbation induced by an external operator 23 by 
considering separately two effects. The first one is the action of the external operator on the unpaired nucleon, while 
the doubly-magic core remains unperturbed. The second effect considers the interaction of the external probe with 
one of the nucleons of the doubly magic core. In this case, the whole nucleus responds to the external perturbation 
because the nucleon struck by the external probe interacts strongly with all the other nucleons, included the unpaired 
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one. 

We express the expectation value of the external operator 23 between two states of the odd-even nucleus with one 
nucleon more than the doubly-magic one as [4-7] 

{A + a\‘B \ A -\-1] (3) = ($0 I Oa ® I 4^0) + ('ho j Kes ^ A{^ap) 25 j 'I'o), (1) 

where and a are the usual creation and annihilation operators, Vres is the residual nucleon-nucleon interaction, 
and Va is a propagator which describes the excitation of the doubly magic core and depends on — ep, the 

difference between the energies of the two states of the A + l system. A similar equation can be written for the A — 1 
nucleus by exchanging creation and annihilation operators. 

In our approach, the excitation of the doubly magic core is described by using the RPA theory. In terms of an 
excitation operator Q'^, we define the excited states of the A-nucleon system as 

I'h,) = Q+jvho), (2) 


while the relation 


Q. I ^-0) = 0 


( 3 ) 


defines the ground state. The label ly indicates all the quantum numbers needed to characterize the excited state. 
The RPA excitation operator is defined as 


Qi ~ ^ ^ Yph ) 

ph 


( 4 ) 


where the index p refers to s.p. states above the Fermi surface (particles) and h to s.p. states below it (holes). The 
RPA propagator can be written as [4, 7, 15] 




^a/3 


i' “r 


( 4-0 I 0.^, 


( 5 ) 


where uii, is the excitation energy of the state | 4'i,). By using this propagator, we express the second term in the 
r.h.s. of Eq. (1) as 


(^o|aaa+EresT5^^(ea;3)'B|To) = 

= ( ^'o I Oa Ofl Eres ^ <5+ | ^'o ) ( - - 

. V - 

~ 'y { 4^0 I 0,a O'p I^es j Qti I 4*0 ) 


1 


~ Wjy Cap -f U!i> 

I I 


^olQ.Sj^o) (6) 

{^o\ [Q.,23] |4/o). 


where we have exploited Eq. (3) to insert the commutator [, ] between operators, with the aim of using the quasi¬ 
boson-approximation [16]: 


{^o\ac,a+VresyT^{ec.p)'B\^o) - 

— ^ ^ ( 4^0 I ^p I^es 1 Qi/ I 4^0 ) 


( 7 ) 


I 


I 




(Choi [Q.,®] jcho). 


The operator describing the action of the external probe on the nuclear system can be expressed in terms of a 
multipole expansion. Each term of this expansion, ®jm, is characterized by the angular momentum J, and its 
projection M on the quantization axis. In our calculations we consider the spherical symmetry of the problem. We 
use a set of s.p. wave functions with the following angular momentum coupling structure: 


<(’o!(r) - ^ ^ Mo! „ Sfl j Jo! Wq ) (d, (ji) Xs 


( 8 ) 


where (r, 0, (f)) are the usual polar coordinates, Yi^ are the spherical harmonics, x the Pauli spinor, and the symbol (|) 
indicates the Clebsch-Gordan coefficient. We calculate the matrix elements by applying the Wigner-Eckart theorem 
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with the phase conventions of Ref. [17], and we sum on all the possible values of rua and mp. For the matrix element 
in Eq. (1) we obtain 


(A + 1; a II 25j 11 A + 1;/3) — 

= {<i^a\\'S>j\\(j)p) + 5a,P (</>/» II 23j II ) 


(9) 


h 

Y.[v''apY.^ph(^p\\'^j\\M + S-^apY. (-.a;];, ( 11 ® J 11 </)„ ) 

ph 


ph 


where the two kernels are defined as 


^ ^ ^a^ph 
ph 


X, 


ph 


^0.8 ^1/ 




f-aP + Wi, 


and 


sss = i: 

ph 


viphp (+ (-1^+^” 


^ph 


€q:/3 + UJy 


^ol 8 


( 10 ) 


( 11 ) 


with 


■^ph = x;^ + . ( 12 ) 

In Eq. (9) the double bars indicate that in the matrix elements between the two s.p. wave functions the angular part 
is evaluated in terms of reduced matrix elements. The calculation of the radial integrals is understood. 

The expression of the interaction terms is analogous to that of the RPA [15] 


^ap'yS 


Y (-1)^+1'’+^^ K 

K 


ja jp J 1 
js jj K J 


X [( Ja K II R,es II 3P 35 K) - (-iyP+^^-^ ( Jg Jp K II Kes || Ja 3^ K ) 


(13) 


In the above equation we used the Wigner 6-j symbol [17], and j = ^2j + 1, where j is referred to an angular 
momentum. The interaction Id-es is described as 




Vihij) + V2{nj)T{i) ■ T{j) + V3{nj)cr{i) ■ (T{j) + V4,{r^j)(T{i) ■ (T{j)T{i) ■ T{j) 

+ tpy-T{i)-T{j)] + 2*Wo X(5(rij)^ij] ■ S + vcouiir^j), 


(14) 


where tp and Wq are real constants, v{rij) are scalar functions of the distance between the two nucleons, and r’coui(Pii) 
is the Coulomb interaction acting between two protons. In the above equations cr and r are the usual Pauli matrices 
for the spin and isospin operators respectively, p is the nucleon density, S = [(t(i) + cr{j)]/2 and p, is the relative 
momentum of the interacting pair with the arrows indicating the side on which the operator acts. It is worth pointing 
out that the same effective nucleon-nucleon interaction is used in the HE and in all the RPA calculations. 

The sum over v in Eq. (9) runs over all the excited states of the A-nucleon core with the same multipolarity and 
parity of the transition operator. Eq. (9) is the basic expression of our model and is analogous to those obtained in 
Refs. [5-7] within the FES theory. 

The results obtained by using the whole expression in Eq. (9) have been labelled as RPA. We compare them with 
those of the IPM, i.e. those provided by the hrst term of Eq. (9). In addition we consider also the results obtained 
within the framework of the traditional time-independent first order perturbation theory (EOPT), which we obtain 
by calculating 


(A-I-1; a II 23 J II A-I-1; /3)fopt = ( 11 ^ j 11 0/3 ) + 5a^p 'Y, ( 11 ^ J 11 <)’^ ) 


E 

ph 


7 / 

^OL^hp 


1 


1 


+ 

^ph ^OL0 H” ^ph 


(0P II II 0/t) • 


(15) 


Eor nuclei with one nucleon less than the doubly magic ones we obtain equations similar to (9) and (15) with an 
additional overall phase. 
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The first operator we consider in this paper is a one-body electric operator, 

-4 


-JM 


-E 


- — ti 
2 




where e is the elementary charge. We express the reduced s.p. matrix element of this operator as 


((/>a II £j||(/'/3) = 


•\/4^ \2 


1 

- -tr 




Y Ip J) ja jp J 


J jp 


-- 0 - 

2 2 


aJ 

‘^a/3 • 


(16) 


(17) 


In the above expression we used the Wigner 3-j symbol [17], and ^(Z) = 1 if Z is even and zero otherwise. We have 
indicated with the integral 


Specifically, we calculate the electric quadrupole moment Q defined as 


n - , f 2 jo 

-Jc 0 Jo 


(A±l;a|| £2||A±1 ;q!). 


(18) 


(19) 


The second operator we consider in this work is the magnetic operator, composed by one- and two-body terms. We 
express the one-body (OB) part of the magnetic operator as 






+ gs{U)(T{i) 


■ V, [r/rjM(0,)] , 


( 20 ) 


where gjsi is the nuclear magneton, and we have for the gyroscopic constants the values gi=l and 0, and g^ = 2.793 
and —1.913, for protons and neutrons, respectively. The expression for the s.p. matrix element is 

{cjo^ms^Up) = ^St^t^i-iy^-^/^+^aic. + ip + j+i)3Jpj (^^1 i) 


X (t/ f^a^) 


gi{to) (l + y^) - 9s{ta) 




( 21 ) 


where we have defined 


— iJoL Ja) (2ja + l) + (l^ J^) (2j^ + 1) . 


( 22 ) 
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A B 

Figure 1; Schematic representation of the MEC diagrams considered in this work. The wavy lines represent the external 
electromagnetic probe, the dashed lines the exchanged charged pion, and the full lines the nucleons. The A diagram is normally 
called seagull and the B diagram pionic. 


In addition to the one-body operator, we consider two-body terms generated by the exchange of a single, charged 
pion, the so-called seagull, M®, represented by the A diagram of Fig. 1, and pionic, represented by the B diagram 
of the same figure. We give in Appendix A the expressions of the matrix elements of these two MEC operators. In 
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the FFS theory, these MFC corrections to the OB operator are taken into account by considering the effective charge 
of the quasi-particle and by including an effective tensor term. 

Specifically, we calculate the magnetic dipole moment defined as 


M = 


f ja 1 ja 
V 3 y ja 0 ja 


(A ± 1; a II Ml II A ± 1; a) , 


(23) 


where Mi = M®® -I- M^®®. In IPM calculations, by considering Mi = we obtain the well known Schmidt 

values [16, 18] which depend only on the angular momentum j of the s.p. level of the unpaired nucleon and are given 
by: 


AiSch 



gi + 9s 


9n , 


for j = / -b i , 


j 

j + 1 



9N , for j 


(24) 


A remarkable difference between our calculations and those based on the FFS theory is that we use bare, not 
effective, electromagnetic operators, therefore we do not include free parameters in our theory. Furthermore, we 
do not include in the magnetic operator (20) the tensor terms often added to take care of magnetic dipole Al = 2 
transitions in the IPM picture [6, 8, 19]. 


III. SPECIFIC APPLICATIONS 

As pointed out in the Introduction, we carried out our calculations by using two different parameterizations of the 
finite-range density dependent Gogny interaction. These define the scalar functions v{rij) of Eq. (15), and the values 
of the tp and Wq constants. We have chosen the traditional DIS force [20] and the more recent DIM one [21], which 
improves the behavior of the neutron matter equation of state at high density values. We would like to emphasize 
the fact that the parameters of these interactions have been chosen once forever in a fit process, described in detail 
in Ref. [13], and they have not been modified in the present work. Since the effective nucleon-nucleon interaction is 
the only input of our calculations, we may state that they are parameter free. 

For each nucleus considered, the calculations of Q and fi are carried out in three steps. In the first one we generate 
the set of s.p. wave functions by means of a HF calculation. In the second step, for a selected multipole excitation 
of angular momentum J and parity 11, we solve the RPA equations. In the third step we calculate the expectation 
value of the operators £2 and Mi by using the expression (9) or the ones for IPM and FOPT. 

From the computational point of view, it is necessary to ensure the numerical convergence of each of these three 
steps. In our calculations we use a set of s.p. wave functions with bound boundary conditions at the edge of an 
r-space box. In this manner all the s.p. states are bound, even those with positive energy. In the common jargon 
these are called discrete calculations. For the first two steps, the numerical stability of the results is related to the 
dimensions of the r-space integration box and to the size of the s.p. configuration space used in the discrete RPA 
calculations. We handle these problems by using the strategy described in Ref. [22], i.e. by choosing the size of the 
integration boxes and the dimensions of the s.p. configuration spaces such as the centroid energies of the electric 
giant dipole resonances do not change by more than 0.5 MeV if the value of any of the two parameters is increased. 
In this respect, the most demanding calculations are those carried out for the ^°®Pb nucleus where we have used a 
box radius of 25 fm and a maximum s.p. energy of 100 MeV. 

The stability of the results of the third step is related to the maximum value of the RPA excitation energy, Amax, 
which is the upper limit of the sum on v in Eq. (9) and of the sum on ph in Eq. (15). We carried out convergence 
tests for each nucleus considered. An example of these tests is presented in Fig. 2, where we show the values of the 
quantities \Q/e\/ZR^, with R = 1.2 in the upper panel, and g/gsch, in the lower one, for the ground states of 
various nuclei, as a function of We observe that in Fig. 2 all the results have already converged at Emax = 70 

MeV, and this is the minimum value of Emax we adopted in our calculations. 


IV. RESULTS 

We have applied the model presented in the previous sections to various nuclei in different regions of the nuclear 
chart. The core nuclei have been selected to have closed shells and spherical symmetry, therefore, in these nuclei, 
deformation and pairing effects do not play any role [23, 24]. We have identified isotopes of oxygen, calcium, zirconium. 
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Figure 2: Convergence test of the stability of our results against -Emax, the upper limits of the sum on v in Eq. (9), for Q, 
panel (a), and jj., panel (b), as given by Eqs. (19) and (23), respectively. The vertical lines indicate the 70 MeV values which 
we have adopted as the value of Emax in our calculations. 


tin and lead endowed with these characteristics. We have investigated the role of the pairing in these nuclei by carrying 
out calculations in a HF plus Bardeen Cooper and Schrieffer framework [24]. The protons and neutrons energy gaps 
between particle and hole states are so large that we did not find pairing effects. We found only one exception, the case 
of the protons in ®°Zr, where the particle fluctuation number is not exactly zero. The evaluation of the occupation 
probabilities of the single particle levels indicates deviations of few percent from that predicted by a sharp Fermi 
distribution. Only the proton 2pi/2 state, which in IPM does not contributes to the Q, is sensitively affected by the 
pairing with an occupation probability of 85%. This s.p. state does not contributes to the Q. 

We first present the results obtained for Q, and, in a second step, those related to fx. 


A. The 2'^ excitations and the Q values. 

We show in Table I the excitation energies and the B(E2)t values of the first 2+ state obtained in our RPA 
calculations for the various nuclei studied, and we compare them with the available experimental values taken from 
the compilations of Refs. [25-27]. We notice a large difference between the RPA values and the experimental ones 
in and '^*^Ca nuclei. The experimental energies are smaller than those predicted by our calculations, and the 
experimental B(E2)t values much larger, especially in ^°Ca. Evidently, for and “^^Ca, the lowest 2+ states cannot 
be properly described by the RPA approach since their structure is more complex than a simple combination of 
one-particle one-hole excitations. In effect, the experimental excitation spectrum of presents at least other six 2+ 
states up to 15 MeV, and that of "‘^Ca about fifty of them below 10 MeV [25, 26]. 

The situation for the other nuclei is different. We found reasonable agreement with the available empirical data 
in 48Ca, i32Sn and ^ospb. The 2+ states in these nuclei are dominated by particle-hole (ph) transitions between 
bound s.p. states. The simplest case is that of the ^®Ca nucleus where we found that the neutron (2p3/2l/^2) 
transition is the dominant component of the first 2+ excited state. We identify analogous situations in ^^^Sn, with 
the proton { 2 d^/ 2 ^gQ/\) and neutron ( 2 / 7/2 l^]"i/ 2 ) transitions, and in ^°®Pb with the proton ( 2 / 7/2 and 

( 1 ^ 9/2 l^]”i/ 2 ) and neutron ( 259 / 21 ^/ 3 / 2 ) transitions. In these nuclei, the 2'^ strength is concentrated in the s.p. 
transitions above mentioned with little fragmentation produced by the coupling with the continuum and/or with 
more complex excitations modes. 

The role of the continuum in this type of calculations has been studied in Ref. [28] where a detailed discussion 
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Table I: Excitation energies, uj and B(E2)'|' valnes of the lowest 2^ state for the core nuclei considered. The results obtained 
with the two interactions we have adopted are shown. The experimental valnes are taken from Refs. [25-27]. 


about the position of the centroid energies and the sum rule exhaustion of the 2'^ strength distribution is presented, 
resulting compatible with the theoretical expectations [29]. There is a difference between the results of Ref. [28] and 
those of the present work, since in these latter RPA calculations the spin-orbit and Coulomb terms of the effective 
nucleon-nucleon interaction have been considered by using the procedure of Ref. [14]. We have verified that for the Q 
values, and also for those of ji, the differences between the two types of calculations are of the order of few percents. 

If we neglect the cases of the and ^°Ca nuclei, the average relative difference between our excitation energies 
and the experimental ones is about 15%. Even though our RPA description of the low-lying 2+ excitation of the core 
nuclei is not satisfactory, we have to consider that, in the calculation of Q in the odd-even nearby nuclei, the relevant 
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^^11/2 

25.62 

27.08 


"Sc“ 

1/5/2 

-10.80 

-10.36 


207 

2^572 

15.59 

16.45 


39Ca 

^^3/2 

2.20 

2.23 


209Bi 

1^9/2 

-27.71 

-29.11 


"Ca 

1/7/2 

-3.68 

-3.66 


209Bi“ 

2/7/2 

-21.31 

-22.28 


47k 

^^3/2 

5.90 

5.97 


209b;6 

1*13/2 

-29.27 

-30.79 


"K“ 

l‘^5/2 

7.34 

7.43 


207pb 

2/5-/1 

11.62 

11.40 



1/7/2 

-10.32 

-10.45 


207pb“ 

^P3/\ 

7.49 

7.37 


«Sc“ 

1/5/2 

-10.52 

-10.63 


207pbb 

^U3/2 

18.54 

18.58 


"Ca 

l/,-/2 

4.95 

4.86 


207pb<= 

2 / 7/1 

12.97 

12.80 


^9Ca 

2p3/2 

-2.58 

-2.40 


207pbd 

^^9/2 

19.53 

19.80 

O 

o 

59k 

^*^3/2 

4.28 

4.61 


209pb 

^99/2 

-14.42 

-13.93 


=9k“ 

11^5/2 

5.67 

6.08 


207pb 

2/5/2 

11.62 

11.40 


6iSc 

1/7/2 

-8.17 

-8.70 


2«9pb“ 

lbl/2 

-21.80 

-21.58 


61Sc“ 

1/5/2 

-7.63 

-8.16 


209pbb 

lil5/2 

-20.11 

-19.81 

9“Zr 

®iNb 

199/2 

-15.55 

-15.94 









8.04 

8.01 







91 Zr 

I97/2 

-8.62 

-8.48 







Table II: Values of Qrpa, in efm^, obtained in our model for the two interactions considered. We present the core nuclei, the 
associated Ail nuclei and the s.p. states which characterize them. The upper latin letters indicate the excited states of the 
isotope. 
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quantity to be considered is the global strength distribution. Our RPA calculations generate the expected amount of 
2+ strength, but it is not properly distributed. This may affect our Q results because the energy denominators of the 
kernels (10) and (11) indicate that the low-lying excited states are more important than those with high energies in 
the evaluation of the matrix element in Eq. (9). We have tested this sensitivity by artificially changing of the 10% 
the values of the uii, energies calculated with the DIM force. The largest deviations of the Q values we have identified 
are of about 10%, much smaller in any case than the differences with the 1PM values. 

Around the core nuclei, listed in Table 1, we have selected 48 states of several isotopes and isotones with one nucleon 
more or less, which allow us to compare our results with the experimental data of the compilation of Ref. [30] and 
also with the results of Ref. [5] obtained within the FFS theory. 

In Table 11, we show the RPA results obtained in our model by using the DIM and DIS nucleon-nucleon effective 
interactions. In this table we indicate the core nuclei, the specific isotopes for which Q was calculated, and the s.p. 
states characterizing them. In general, we have considered the states with lowest energy and Q ^ 0. For example, 
the true ground states of ^°^Pb are a 1/2+ and a 1/2“ states respectively, but both have Q = 0. For this 

reason, in these two cases, we have indicated the 3/2+ and 5/2“ states, respectively. In the table, the states of the 
same isotope with higher energy are indicated by a superscript latin letter. 

Our calculations confirm the well known fact that the Q values of nuclei with one nucleon more than the doubly 
closed shell core are negative, and those with a nucleon less are positive [19]. A first general remark about the results 
of Table 11, is that the values we have obtained by using the two different interactions are very similar. The average 
relative difference between the RPA results calculated with DIM and DIS interactions is of few percent. Henceforth, 
if not explicitly stated, we shall discuss only the results obtained with the DIM interaction. 

We compare in Fig. 3 the results of the 1PM, triangles, FOPT, squares, and RPA, circles, calculations. The four 
panels show separately the results obtained for hole and particle-like states and for protons and neutrons. Since we 
do not use effective charge, the 1PM results for the neutron states are zero. 


protons 


0 “ 



.».IPM 

23 41 49 61 ^91 B3 209”^ 

17 25 41“ 49“ 61“ 101 209 209* 



0 

'Ca 

Zr 

'Sn' 

Pb ( 

C) 

- ^ 





- 











•m -- 


particle states 



- 


41 

9i 

133 



77 

49 
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209 209 



Ca 

Zr 

Sn 

Pb 




hole states 





- Z'"A' 


'A'"A" 


47 


99 ^ 4 

07 207*^207'' 

39 

89 

131 

207" 207^ 


Figure 3: Values of Q, in e fm^, for RPA (circles), FOPT (squares), and IPM (triangles), calculated with the DIM interaction 
for the various nuclei studied. The left (right) panels correspond to proton (neutron) states, while the upper (lower) panels 
show the results for particle (hole) states. The lines have been drawn to guide the eyes. The latin upper indexes indicate 
excited states following the nomenclature of Table II. The RPA values for the nuclei around the ^^®Sn core, drawn as open 
circles, have been multiplied by 0.25. 


A first, general, remark is that, in absolute value, the RPA results are always larger than those obtained in the 
other calculations. It is interesting to observe that, in the proton case, the FOPT results are very close to those of the 
IPM, and, sometime, even smaller in absolute value. In general, the RPA results produce a relatively small correction 
with respect to the IPM values. 

We carried out an additional test of our model by comparing our results with those of other theoretical approaches, 
the Hartree-Fock-Bogolioubov (HFB) and the FFS theory [7, 10-12]. This comparison is shown in Fig. 4. The HFB 
calculations have been performed by using the technique developed by Robledo et al. [31], based on the gradient 
method. In the minimization procedure both time-even and time-odd fields are considered in the calculation of the 
energy functional and this allows us to avoid the equal filling approximation [32]. In the figure, we show the RPA and 
HFB results obtained with the DIS interaction. The relevant point for the present study is the remarkable agreement 
between the HFB and our RPA results. 

The other comparison proposed in Fig. 4 is that with the results of FFS calculations. This last approach is 
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Figure 4: Values of Q, in e fm^, obtained in our RPA calculations with DIS interaction (solid circles), compared to the FFS 
results taken from Refs. [7, 10-12] (open squares) and to the HFB results found for DIS interaction (open circles). 


rather similar to ours, and the differences are due to the use of different interactions, and s.p. bases. The results of 
the nuclei around taken from Ref. [7], have been obtained by fully exploiting the philosophy of the Landau- 

Migdal approach. The s.p. wave functions have been generated by using a Woods-Saxon potential well, and the 
RPA calculations have been conducted with a Landau-Migdal interaction whose parameters have been chosen to 
reproduce at best the data. The results of nuclei lighter than ^°®Pb, are taken from Refs. [10-12], and have been 
obtained by using a slightly different computational scheme. Skyrme interactions of different type have been used in 
HF calculations to generate the s.p. basis wave functions and, in the second step, also to solve the RPA equations. 
The zero-range characteristics of this type of interaction requires the use of a cutoff renormalisation parameter whose 
value depends on the size of the s.p. configuration space used in the RPA calculations. The Q values obtained with 
the FFS approach are, in absolute value, larger than ours. 

In Table III we compare our results with the experimental values taken from the compilation of Ref. [30] and with 
those given in Ref. [33]. The RPA description of these experimental data is not particularly satisfactory. All the 
experimental values are larger than our predictions. In heavier nuclei our calculations are able to account only for 
about half of the observed Q values. This indicates that our approach is not fully able to describe the large collectivity 
of the 2+ excitations. On the other hand, the improvement with respect to the IPM predictions, and also with respect 
to the FOPT predictions, is remarkable. 


exp 


core 

nucleus 

s.p. state 

IPM FOPT 

RPA Ref. [30] Ref. [33] 

160 

ITp 

1^5/2 

-4.12 

-3.74 

-4.50 

-5.80 



i"o 

1^5/2 

0.0 

-1.51 

-1.61 

-2.60 


o 

o 

39k 

^^3/2 

3.37 

3.04 

4.33 

5.90 



«Sc 

1 / 7/2 

-7.01 

-6.31 

-8.30 

-15.60 



3®Ca 

1^3/2 

0.0 

1.78 

2.20 

3.60 



«Ca 

1 / 7/2 

0.0 

-3.17 

-3.68 

-9.00 


«Ca 

47Ca 

1/772 

0.0 

2.43 

4.95 

2.10 

8.46 


«Ca 

2P3/2 

0.0 

-1.24 

-2.58 


-3.63 

®0Zr 

Zr 

2d5/2 

0.0 

-3.15 

-5.68 

-17.60 


208 Pb 

209Bi 

1^9/2 

-16.39 

-13.84 

-27.71 

-51.60 



209pb 

2(?9/2 

0.0 

-7.55 

-14.42 

-30.00 



Table III: Values of Q, in efm^, for IPM, FOPT and RPA, calculated with the DIM interaction, compared to the experimental 
data taken from the compilation [30] and from Ref. [33]. The s.p. states refer to the corresponding doubly magic cores. 
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B. The fi and the 1’^' excitation. 

We use the same set of even-even core nuclei adopted for the calculation of Q also for the investigation of the 
magnetic dipole moment /i. Around each core nucleus we calculate the /r values for the four nuclei with one nucleon 
more or less. We present our results by following the strategy adopted in the investigation of Q. We first discuss 
the features of the I"*" excitation in the core nuclei and, after, we present the results of the calculation of ^ and we 
compare them with the experimental values. 

We present in Table IV the excitation energies and the B(Ml)t values of the main low-lying 1+ states, and we 
compare them with the experimental values taken from Refs. [25, 26]. We indicate in the table also the dominant ph 
transitions. 


nucleus 

UJ 

(MeV) 

B(Ml)t (Mi) 

main s.p. 

transitions 

DIM 

DIS 

exp 

DIM 

DIS 

proton 

neutron 

16 o 

17.76 

18.37 

13.66 

0.01 

0.01 



22 0 

8.35 

8.97 


5.81 

5.18 


(l(l3/2l(l5/2) 

240 

8.50 

9.06 

9.50 

5.75 

5.13 


(ld3/2l(l5/2) 

40Ca 

14.62 

15.03 

9.87 

0.006 

0.006 

(2^5/21(13/2) 


48 ca 

9.30 

10.19 

10.23 

10.23 

9.66 


(1/5/2 1/7/2) 

80 Ca 

6.73 

6.76 


0.04 

0.04 


(3p3/22p“ 2 ) 

90 Zr 

9.08 

9.98 

9.37 

13.98 

13.49 


(I57/2I99/2) 

lOOSn 

7.49 

9.13 


0.65 

0.97 

(l/ 5 / 2 l/ 7 “/ 2 ) 


l32Sn 

6.78 

8.00 


6.05 

9.84 



208 pb 

6.30 

7.60 

5.85 

7.80 

11.89 

( 1 ^ 11 / 21 ^ 13 / 2 ) ( 1 ^ 9 / 21 ^ 11 / 2 ) 


Table IV: Energies, cu, and B(Ml)b values of the main low-lying 1+ excited state obtained in our calculations for the core nuclei 
considered, with the two interactions we have adopted. The experimental values are taken from Refs. [25, 26]. 

The table shows that the B(Ml)t values of the ^®0, ^°Ca and ®°Ca nuclei, where all the spin-orbit partner levels 
are occupied, are orders of magnitudes smaller than the other ones. This confirms the well known fact [19, 34] that 
the Ml excitation is strongly excited in nuclei where ph transitions between spin-orbit partner levels are allowed. We 
have presented in Ref. [35] a detailed discussion of the Ml strength distribution in oxygen and calcium isotopes. 

The agreement with the experimental energies is reasonable for almost all the cases, but for and ^°Ca. We have 
already pointed out the fact that in the first two nuclei the I"*" excitation is rather weak because all the spin-orbit 
partner levels are occupied. In these nuclei it is not easy to find out the I’*' states to be compared with those predicted 
by our model, since, experimentally, the Ml strength is rather fragmented. For example in ^'^Ca we found more than 
twenty 1+ states below 12 MeV [25, 26]. 

Our description of the 1^ excitation is rather good in the situations where these excitations are dominated by a 
few ph transitions. In these cases our approach works at its best. The main features of the transitions are already 
well described within the IBM which is corrected by the RPA by including other, less important, ph transitions. 

The ten doubly magic nuclei listed in Table IV have been considered as core nuclei where we add or remove one 
nucleon at the time to form odd-even nuclei with magnetic dipole moment /r. Before discussing the effects of the core 
polarization, i.e. of the excitation of the even-even core, we clarify the role played by the MFC in the calculation of 

fi. 

In Fig. 5 we adopt the usual representation of /r [16] against the angular momentum j of the odd-even nucleus, 
to show our results obtained with the DIM force. In the (a) and (b) panels the symbols represent the /iipM values 
obtained by considering the full electromagnetic operator Mi = M®^ -I- while the lines indicate the Schmidt 

values found for M®® (see Eq. (24)). The MEG produce relatively small changes with respect to the Schmidt 
values, usually of the order of few percent up to a maximum value of about 10%. These changes have a consistent 
behavior over all the various regions of the nuclear chart, and we did not observe remarkable differences between 
nuclei characterized by a hole- (open squares) or a particle-like (solid circles) structure. 

Another perspective to observe these results is given in Fig. 6 where we show the differences ^ipm — MSch as a 
function of the mass number of the nuclei investigated. The results of panel (a) correspond to the odd-proton nuclei 
while those of panel (b) are those obtained for the odd-neutron nuclei. In the odd-proton cases the MEG always 
increase the /r values with respect to the Schmidt ones, but for the and nuclei. The results of panel (b) 
show that the effect of the MEG on the odd-neutron nuclei consists in a lowering with respect to the Schmidt values, 
without exceptions. It is possible to outline a general overall tendency of the MEG effects to increase with the mass 
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IPM 



RPA 




13 5 7 

j + 1/2 


Figure 5: Comparison between the values of /tipm, left panels, and /trpa, right panels, expressed in nuclear magnetons, and 
the Schmidt values, as a function of the angular momentum j of the odd-even nucleus. The results obtained with the full 
electromagnetic operator. Mi = M?® -I- are indicated by the symbols and the Schmidt values by the solid lines. Open 

squares and solid circles represent the results obtained, respectively, for nuclei characterized by hole or particle-like s.p. states. 
The calculations have been carried out with the DIM interaction. 


number, even though there are quite a few exceptions. 

The effects we have just described are consistent for all the nuclei, and, in terms of the agreement with the 
experimental data implies two contrasting aspects. Since the experimental data are lying between the two Schmidt 
lines, the inclusion of the MEG improves the agreement for odd-proton nuclei with j< = I — 1/2, while they are 
worsening it for odd-proton nuclei with = I + 1/2. We have a similar situation also for the odd-neutron nuclei 
where we observe an improvement in the case of j< and a worsening for />. 

In Table V we present the RPA results obtained with the two interactions we have adopted, and by using the 
full electromagnetic operator which includes the MEG. We have considered those configurations which generate the 




Figure 6: Difference pipM — /rsch, expressed in nuclear magnetons, for proton, panel (a), and neutron, panel (b), odd nuclei. 
Open squares and solid circles represent the results obtained for nuclei characterized by hole or particle-like structure. The 
calculations have been carried out with the DIM interaction. 
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core 

nucleus 

S.p. state 

DIM 

DIS 

core 

nucleus 

S.p. state 

DIM 

DIS 

16o 

15n 

^Pl/2 

-0.137 

-0.142 

“Zr 

89y 

207/2 

-0.108 

-0.152 


i^F 

ld5/2 

4.884 

4.885 



lg9/2 

6.727 

6.996 


ISq 

ipr/2 

0.638 

0.638 



^89/2 

-1.937 

-1.813 


1^0 

1^5/2 

-2.088 

-2.007 


"Zr 

26-5/2 

-1.808 

-1.715 

220 

21n 

IPfA 

0.001 

-0.042 

i™Sn 


^89/2 

6.559 

6.593 


23p 

1^5/2 

4.397 

4.779 


i"Sb 

265/2 

4.448 

4.476 


21o 

1^7/2 

-1.667 

-1.487 


'>'>Sn 

^89/2 

-1.638 

-1.771 


23o 

2Si/2 

-1.610 

-1.485 


i"Sn 

265/2 

-1.507 

-1.614 

240 

23n 

ipr/2 

0.013 

-0.024 


i"ln 

^89/2 

6.788 

6.767 


25f 

1^5/2 

4.444 

4.801 


i33Sb 

lg7/2 

2.579 

2.546 


23o 

2s^ ^ 

-1.628 

-1.482 


i"Sn 

^'^3/2 

0.684 

0.758 


25o 

1^3/2 

0.842 

0.841 


l3lsn- 

^^11/2 

-1.520 

- 1.688 

^“Ca 

39k 

^'^3/2 

0.358 

0.350 


i^^Sn 

169/2 

0.596 

0.746 


"Sc 

1^7/2 

5.971 

5.972 

208pb 

207 rpj 

1/2 

2.528 

2.538 


3^Ca 

^^3/2 

0.910 

0.918 


209Bi 

169/2 

3.604 

3.548 


"Ca 

1^7/2 

-2.096 

-2.098 


2“9Bi“ 

lil3/2 

8.859 



47k 

2 s~^ 

•^^1/2 

2.384 

2.695 


207pb 

3p7/2 

0.437 

0.463 



lf7/2 

5.577 

5.888 


2“7pb“ 

^*^5/2 

0.764 

0.857 


"Ca 

lf-1 

-^^7/2 

-1.809 

-1.664 


207pbb 

^' 13/2 

-1.655 



^s>Ca 

2p3/2 

-1.709 

-1.605 


209pb 

2g9/2 

-1.635 

-1.756 

8“Ca 

59k 

^■^3/2 

0.347 

0.338 







"Sc 

1^7/2 

6.146 

6.147 







S^Ca 

^P3/2 

-1.954 

-1.951 







"Ca 

lg9/2 

-2.015 

- 2.011 







Table V: Values of fi, in nuclear magnetons, obtained in our RPA model with the full operator for the two 

interactions considered. We present the core nuclei, the specihc Ail nuclei and the s.p. states which characterise them. The 
upper latin letters indicate the excited states of the isotope. 


ground state of the odd-even nuclei, and some excited state of interest. In the table, our results are compared to the 
available data taken from the compilation of Ref. [30]. 

With the exceptions of the and ^^^Sn cases, some of which we shall discuss in more detail, the relative 

differences between the DIM and DIS results are of the order of 1%. This indicates that the effects we are discussing 
are almost independent of the specific implementation of the interaction used. 

The quantitative stability of our results can be observed by considering that, in our approach, the nucleus 
can be obtained by adding a neutron to the core or subtracting one to The differences between the two 
procedures are related to the different s.p. states produced by the HF calculations and the results obtained for DIM 
and DIS interactions show /i values that differ on the third and fourth significant figures, respectively. 

We compare the DIM results of Table V with the Schmidt values in the panels (c) and (d) of Fig. 5. In general, the 
Z^RPA values are situated between the two Schmidt lines, as it is observed for the experimental values. The deviation 
from the Schmidt values is particularly evident in the neutron case for the nuclei with j< above 5/2, which correspond 
to the i33Sn isotopes. 

In Fig. 7 we compare the RPA, FOPT and IPM results. All the calculations have been carried out with the DIM 
interaction, and by using the complete electromagnetic operator Mi = Mf® -I- M^^^. The results are shown as 
relative difference with respect to the Schmidt values, Eq. (24). We observe the great similarity between our RPA 
results and those obtained in FOPT. This indicates that the structure of the 1+ states is not as collective as that of 
the 2+ states, but it is rather dominated by single-particle excitations. 

The FOPT and RPA results show very similar behaviors. In general, the FOPT deviate more from the Schmidt 
values than the RPA ones. We remark that the FOPT approach, often utilized in the literature [3, 36-39], performs 
much better in this case than for the electric quadrupole moment case, as it is shown in Fig. 4. 

Our results agree with those of the non-relativistic calculations of Refs. [3, 38]. More precisely, and always with 
respect to the Schmidt values, we observe that: (i) both MEC and RPA increase the /i values in ^°®Bi; (ii) in 207qq .^g 
find an enhancement due to the MEC while the RPA effects produce a lowering of /i; (iii) in ^°®Pb we have opposite 
effects, and (iv) in ^°^Pb both MEC and RPA lower the /i values. 

The fi values obtained in RPA show a remarkable difference with respect to the IPM results. The case of the 
and nuclei is particularly interesting since the three calculations provide very similar results and all of them 
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A 


Figure 7: Comparison between our RPA, circles, FOPT, squares, and IPM, triangles, results. The calculations have been 
carried out with the DIM interaction and by using the full electromagnetic operator Mi = Mi +Mi . The lines have been 
drawn to guide the eyes. 


noticeably deviate from the Schmidt values. This is only due to the presence of MEC in the operator. This effect 
should be investigated better, by considering also that the nucleus does not present this feature. 

A direct comparison between our results and the available experimental data taken from Ref. [30] is presented in 
Table VI. In this table we also include the results of the FFS calculation of Ref. [8]. We observe that the effects 
of the core polarisation become more important in the heavier nuclei. The differences between the results of IPM, 
FOPT and RPA calculations are rather small in the nuclei around and ^°Ca. 

In the great majority of the cases the RPA results are closer to the experimental values than the other ones, 
indicating that the inclusion of the core polarization improves the agreement with the available experimental data. 
However, this is not a universal behavior since there are very specific cases, where the IPM approach produce a better 


core 

nucleus 

s.p. state 

Schmidt 

IPM 

FOPT 

RPA 

FFS 

exp 

16o 

15n 

^Pl/2 

-0.264 

-0.153 

-0.133 

-0.137 

— 

-0.283 


"F 

ld5/2 

4.793 

4.896 

4.882 

4.884 

— 

4.721 


ISq 

lPf /2 

0.638 

0.524 

0.506 

0.510 

— 

0.720 



1^5/2 

-1.913 

-2.018 

-2.003 

-2.088 

— 

-1.894 

"«Ca 

39k 

^'^3/2 

0.124 

0.340 

0.361 

0.358 

0.329 

0.391 


"Sc 

lf7/2 

5.793 

5.983 

5.963 

5.971 

5.485 

5.431 



^‘^3/2 

1.148 

0.928 

0.910 

0.910 

0.888 

1.022 


"Ca 

1^7/2 

-1.913 

-2.109 

-2.088 

-2.096 

-1.626 

-1.595 

«Ca 

47 k 

2^/2 

2.793 

3.011 

2.173 

2.384 

— 

1.933 


"Ca 

lf7/2 

-1.913 

-2.112 

-1.732 

-1.809 

-1.556 

-1.380 


«Ca 

2p3/2 

-1.913 

- 1.954 

-1.654 

-1.709 

-1.156 

— 


89y 

^Pl/2 

-0.264 

-0.221 

-0.114 

-0.108 

-0.125 

-0.137 



^ 89/2 

-1.913 

-2.212 

-1.872 

- 1.937 

-1.304 

-1.046 


"Zr 

2d5/2 

-1.913 

-2.022 

-1.769 

-1.808 

-1.214 

-1.304 



lg7/2 

1.717 

1.987 

2.826 

2.579 

2.693 

3.000 



^*^3/2 

1.148 

1.029 

0.601 

0.684 

0.681 

0.747 



lhii/2 

-1.913 

-2.191 

-1.227 

-1.520 

— 

-1.276 

208 pb 

20T 

3s-' 

2.793 

2.990 

2.379 

2.528 

1.857 

1.876 


209Bi 

1^9/2 

2.624 

3.025 

3.827 

3.604 

3.691 

4.110 


207pb 

3Pl/2 

0.638 

0.549 

0.472 

0.437 

0.473 

0.593 


2«7pb“ 

^^5/2 

1.366 

1.159 

0.677 

0.764 

0.720 

0.800 


209pb 

2g9/2 

-1.913 

-2.110 

-1.486 

-1.635 

-1.337 

-1.474 


Table VI: Comparison between p values, expressed in nuclear magnetons, calcnlated in IPM, FOPT and RPA by using the 
complete electromagnetic operator M°® with the Schmidt values (24), the FFS results of Ref. [8] and the experimental 

data of Ref. [30]. 
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nucleus interaction Schmidt 

IPM 

FOPT 

RPA 

21 N 

DIM 

DIS 

-0.264 

-0.264 

-0.341 

-0.339 

0.030 

0.003 

0.001 

-0.042 

23n 

DIM 

DIS 

-0.264 

-0.264 

-0.323 

-0.322 

0.025 

-0.002 

0.013 

-0.024 


Table VII: Values of /x, in nuclear magnetons, of the and nuclei. The Schmidt values are compared to those found 
within IPM, FOPT and RPA approaches with the full MP® + operator. 


description of the observed /x values. As expected the FFS results, which use effective electromagnetic operators, are 
in general closer to the experimental values that ours. 

We conclude this section by analyzing in detail the cases of the and nuclei. In Table VII we show the /i 
values for these two nuclei calculated by using different approximations. We observe that the Schmidt values have 
negative sign. The inclusion of the MFC in the IPM calculations gives a negative contribution, further lowering these 
values. The core polarization plays an important role as the results of two last columns indicate. In both nuclei, the 
effect of the residual interaction has opposite sign with respect to the contribution of the MFC. This effect is larger 
in the FOPT than in RPA results. In RPA calculations, the effects of the DIM force are large enough to change the 
sign of /X with respect to the Schmidt value. Fxperimental ^ values for these two nuclei are not available. In any case, 
our RPA results follow the trend of the know experimental values of lying between the Schmidt values. 


V. CONCLUSIONS 

In this work, we have presented the results of a parameter free theoretical approach which describes the properties 
of odd-even nuclei with one nucleon more, or less, than a doubly magic one. The convergence of the energy sum in Fq. 
(9), the fundamental equation of our model, is ensured by the use of finite-range nucleon-nucleon effective interaction. 
The use of zero-range interactions, as it is done for example in Refs. [5-7], requires the inclusion of additional free 
parameters. The universality of our approach allows us to apply it in each region of the nuclear chart, from light to 
heavy nuclei. 

We have considered two parameterizations of the finite-range effective Gogny nucleon-nucleon interaction taken from 
the literature, the DIS and DIM. Our calculations are fully self-consistent, meaning that all results have been obtained 
by using the same interaction in each of the three steps of the calculation: the generation of the s.p. configuration 
space by means of a HF procedure, the solution of the RPA equations and, finally, the evaluation of the response of 
the odd-even nucleus to the external probe. In each step of the calculation the complete interaction has been used, 
including the spin-orbit and the Coulomb terms, commonly neglected in RPA. 

We have tested the validity of our approach by calculating the values of the electric quadrupole moment, Q, of 48 
different states, and the magnetic dipole moment, n, of 44 states, in various regions of the nuclear chart. We have used 
the traditional one-body operator for the calculation of Q, while for the ^ we have considered also the contribution 
of MFC. We use bare electromagnetic operators whose coupling constants are those of free nucleons. 

In absolute value, the results obtained by considering the MFC are always larger than the Schmidt values, implying 
an improvement of the agreement with experimental data for odd-proton nuclei with and with odd-neutron nuclei 
with j>, and a worsening for the other type of nuclei. 

The differences between the results obtained with the two forces are small, negligible if compared to the differences 
between RPA and IPM results. This indicates that the quantitative description of the core polarisation effects in our 
calculations is more related to the theoretical approach, rather than to the use of a specific interaction. 

Our calculations require the description of the excitation of 2+ and I'*' states in doubly closed shell nuclei, in the 
framework of the RPA theory. We observe that our model describe better the excitations than the 2+ ones. 
This occurs because the most important I’*' states, those where the largest part of the strength is concentrated, are 
dominated by a single ph transition. In this case, our approach works at its best, by correcting the main ph transition 
with the presence of other one-particle one-hole transitions. On the other hand, the structure of many of the 2+ states 
is very collective, and a good description of it requires to go beyond the linear combination of one-particle one-hole 
transitions which is what the RPA considers. 

This difference in the description of the two multipolarities is reflected in the results we have obtained for the values 
of Q and The description of the magnetic dipole moments is certainly better than that of the electric quadrupole 
moments. The similarity between our results and those obtained in FOPT is a further indication of the perturbative 
structure of the I"*" excitations. In any case, we have to point out that our results show a remarkable improvement of 
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the description of the experimental values with respect to the IPM predictions, and also with respect to the FOPT 
for Q. 

We found a good agreement between our results and those of the FFS theory given in the literature [5-12]. The basic 
theoretical hypotheses of the two approaches are rather similar. More remarkable is the agreement with the Q values 
obtained in the HFB calculations. Our approach and the HFB treat the core polarisation in a completely different 
manner, but they produce similar results. In our approach the IPM response is corrected by the core polarisation 
described in terms of particle-hole excitations. Alternatively the BFB theory corrects the IPM by considering the 
effects of the pairing force. These two alternative visions of effects beyond the IPM picture generate very similar 
quantitative results. It would be interesting to identify specific observables which allow us to disentangle the two 
effects. 

Despite the limitations we have pointed out, we think that our parameter free, self-consistent approach is adequate 
to make predictions of the properties of odd-even nuclei nearby the doubly magic-ones in regions of the nuclear chart 
not yet experimentally explored. 


Appendix A: MEC matrix elements 


In this appendix we give the explicit expressions of the matrix elements of the MEC inserted in Eq. (9) to calculate 
the values of /r. We have considered the contribution of the two diagrams indicated in Eig. 1 and called them seagull 
and pionic. We give here the final results of a calculation which describes these contributions in terms of Feynman’s 
diagrams, makes a non-relativistic reduction, a multipole expansion and the use of s.p. wave functions with the 
angular coupling indicated in Eq. (8). 

We call q the modulus of the momentum that the photon is exchanging with the nucleus, and uj the photon energy. 
In the specific case under study we have uj = tap and q = uj/(he). 

The contribution of the seagull term is given by 


( </>a II AIj \\ (pb) — ^(la + Ip + J + 1) (— 1) 


iya-i/2 


jajbJ ( ja jb J 


yj(jTT) V 5 -5 0 
X [G^(q,w) - G^(q,w)] J drr^ jj(qr)j‘^^l j(r) , 


where and G^ are, respectively, the proton and neutron electromagnetic form factors, and we have defined 


^1,1,j(^) “ ^ 2 [^ab,pSh,n Sab,nSh,p] 


(Al) 


h<eF 


xi^aia + lh + Li + l)Ll Kab‘rah(Li,r)Rbir)Rb(r) 

+ + lh + L2 + l)Ll ^ 7bb{L2, r) Ra{r) Rh(r) 1 , 


where = 0.079 is the pion-nucleon coupling constant, is the pion mass, k is defined in Eq. (22), ^(l) = 1 if Hs 
even and zero otherwise, R(r) is the radial part of the s.p. wave function (8), and T is defined as 

7ab(L,x) = J drr^ { ((Ir (^)| 

r°° „ 2 

X / dkk -VT,(k, Cab) jL(kr)jL(kx). (A2) 

JO 

In the above equations we have indicated with Vt^ the pion propagator 


e) 


I 

— £2 _|_ ^2 • 


(A3) 


The contribution of the pionic term is 


( (pa II Atj \\ (f>b) — ^{Ict + Ip + J + ja jb J 7^7(9, w) 


drr"^ jj{qr)7^^lj(r) , 


(A4) 






17 


where 



(AS) 


is the electromagnetic form factor of the pion, with rup = 776.0 MeV the p-meson mass, and 



/I<ep 

X ^ ^ ^{la + + Ll + 1 ) ^(/f, + + L2 + 1 ) L\ L\ 


(A6) 



Li J L 2 

1 -1 0 


X \/^2(7^2 + 1 ) r) r). 
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